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ABSTRACT 

We analyze the stability of g-modes in variable white dwarfs with hydrogen 
envelopes. All the relevant physical processes take place in the outer layer of hydrogen 
rich material which consists of a radiative layer overlain by a convective envelope. 
The radiative layer contributes to mode damping because its opacity decreases upon 
compression and the amplitude of the Lagrangian pressure perturbation increases 
outward. The convective envelope is the seat of mode excitation because it acts as 
an insulating blanket with respect to the perturbed flux that enters it from below. A 
crucial point is that the convective motions respond to the instantaneous pulsational 
state. Driving exceeds damping by as much as a factor of two provided lot c > 1, where 
uj is the radian frequency of the mode and r c ~ 4r t h with r t h being the thermal time 
constant evaluated at the base of the convective envelope. As a white dwarf cools, its 
convection zone deepens, and modes of lower frequency become overstable. However, 
the deeper convection zone impedes the passage of flux perturbations from the base of 
the convection zone to the photosphere. Thus the photometric variation of a mode with 
constant velocity amplitude decreases. These factors account for the observed trend 
that longer period modes are found in cooler DAVs. Overstable modes have growth 
rates of order 7 ~ l/inrj), where n is the mode's radial order and is the thermal 
time-scale evaluated at the top of the mode's cavity. The growth time, 7 , ranges 
from hours for the longest period observed modes (P ~ 20 minutes) to thousands of 
years for those of shortest period (P ~ 2 minutes). The linear growth time probably 
sets the time-scale for variations of mode amplitude and phase. This is consistent with 
observations showing that longer period modes are more variable than shorter period 
ones. Our investigation confirms many results obtained by Brickhill in his pioneering 
studies of ZZ Cetis. However, it suffers from at least two serious shortcomings. It 
is based on the quasiadiabatic approximation that strictly applies only in the limit 
ujt c 3> 1, and it ignores damping associated with turbulent viscosity in the convection 
zone. We will remove these shortcomings in future papers. 
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1. Introduction 

ZZ Cetis, also called DAVs, are variable white dwarfs with hydrogen atmospheres. Their 
photometric variations are associated with nonradial gravity-modes (g- modes); for the first 
conclusive proof, see Robinson et al. ( |1983| ). These stars have shallow surface convection zones 
overlying stably stratified interiors. As the result of gravitational settling, different elements are 
well separated . With increasing depth, the composition changes from hydrogen to helium, and 
then in most cases to a mixture of carbon and oxygen. From center to surface the luminosity is 
carried first by electron conduction, then by radiative diffusion, and finally by convection. 

Our aim is to describe the mechanism responsible for the overstability of g-modes in ZZ Ceti 
stars. This topic has received attention in the past. Initial calculations of overstable modes were 
presented in Dziembowski & Koester ( [1981|) , Dolez k. Vauclair ( 1981 ), and Winget et al. ( 1982| ). 



These were based on the assumption that the convective flux does not respond to pulsation; this is 
often referred to as the frozen convection hypothesis. Because hydrogen is partially ionized in the 
surface layers of ZZ Ceti stars, these workers attributed mode excitation to the K-mechanism. In so 
doing, they ignored the fact that the thermal time-scale in the layer of partial ionization is many 



orders of magnitude smaller than the periods of the overstable modes. Pesnell (1987) pointed out 
that in calculations such as those just referred to, mode excitation results from the outward decay 
of the perturbed radiative flux at the bottom of the convective envelope. He coined the term 
'convective blocking' for this excitation mechanising Although convective blocking is responsible 
for mode excitation in the above cited references, it does not occur in the convective envelopes of 
ZZ Ceti stars. This is because the dynamic time-scale for convective readjustment (i.e., convective 
turn-over time) in these stars is much shorter than the g-mode periods. Noting this, Brickhill 



( 1983 , 1990 , |1991a , 1991b ) assumed that convection responds instantaneously to the pulsational 
state. He demonstrated that this leads to a new type of mode excitation, which he referred to 
as convective driving. Brickhill went on to presents the first physically consistent calculations of 
mode overstability, mode visibility, and instability strip width. Our investigation supports most 
of his conclusions. Additional support for convective driving is provided by Gautschy et al. ( 1996| ) 



who found overstable modes in calculations in which convection is modeled by hydrodynamic 
simulation. 

In this paper we elucidate the manner in which instantaneous convective adjustment promotes 
mode overstability. We adopt the quasiadiabatic approximation in the radiative interior. We also 



3 This mechanism was described in a general way by Cox & Guili (196£), and explained in more detail by Goldreich 
& Keeley ( |l977| ). 
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ignore the effects of turbulent viscosity in the convection zone. These simplifications enable us to 
keep our investigation analytical, although we appeal to numerically computed stellar models and 
eigenfunctions for guidance. The DA white dwarf models we use are those produced by Bradley 
(1996) for asteroseismology. Fully nonadiabatic results, which require numerical computations, 



will be reported in a subsequent paper. These modify the details, but not the principal conclusions 
arrived at in the present paper. 

The plan of our paper is as follows. The linearized wave equation is derived in §||. In §j3|, 
we evaluate the perturbations associated with a g-mode in different parts of the star. We devote 
§U to the derivation of a simple overstability criterion. Relevant time-scales and the validity 
of the quasiadiabatic approximation are discussed in §|5[ The appendix contains derivations of 
convenient scaling relations for the dispersion relation, the WKB eigenfunction, and the amplitude 
normalization. 



2. Preparations 

We adopt standard notation for the thermodynamic variables pressure, p, density, p, and 
temperature, T. Our specific entropy, s, is dimensionless; we measure it in units of fee/m^, where 
liB is the Boltzmann constant and m p the proton mass. We denote the gravitational acceleration 
by g, the opacity per unit mass by k, the pressure scale height by, H p = p/(gp), and the adiabatic 
sound speed by c s where c 2 = (dp/dp) s . 

A g-mode in a spherical star is characterized by three eigenvalues (n,£,m); n is the number 
of radial nodes in the radial component of the displacement vector, I is the angular degree, and m 
the azimuthal separation parameter .0 G-modes detected in DAV stars have modest radial orders, 
1 < n < 25, and low angular degrees, 1 < £ < 2.0. Their periods fall in the range between 100 s 
and 1200 s. G-mode periods increase with n and decrease with I (see § |A.3| ). These properties 
reflect the nature of the restoring force, namely gravity which opposes departure of surfaces of 
constant density from those of constant gravitational potential. 

Gravity waves propagate where u is smaller than both the radian Lamb (acoustic) frequency, 
Lt, and the radian Brunt- Vaisala (buoyancy) frequency, N. These critical frequencies are defined 
by Lj = l{l + l)(c s /r) 2 , and iV 2 = — g{d In p/ds) p ds/ dr. Electron degeneracy reduces the 
buoyancy frequency in the deep interior of a white dwarf (cf. Fig. 1). Consequently, the uppermost 
radial node of even the lowest order g-mode is confined to the outer few percent of the stellar 
radius. Since all of the relevant physics takes place in the region above this node, it is convenient to 
adopt a plane-parallel approximation with g taken to be constant in both space and time. In place 
of the radial coordinate r, we use the depth z below the photosphere. We neglect gravitational 



4 We ignore m in this paper since DAV's are slow rotators. 
5 The latter is probably an observational selection effect. 
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perturbations; due to the small fraction of the stellar mass that a mode samples, its gravitational 
perturbation is insignificant. A few order of magnitude relations to keep in mind are: p ~ gpz, 
Hp ~ z, dlnp/dz ~ 1/z, and c 2 ~ gz. In the convection zone, N 2 ~ — (v cv /c s ) 2 g/z ~ — (v cv /z) 2 , 
where v cv is the convective velocity. In the upper radiative interior, A^ 2 ~ g/z ~ (c s /z) 2 . 



The linearized equations of mass and momentum conservation, augmented by the linearized, 
adiabatic equation of state, read[] 
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In the above equations, ^ and £ z represent horizontal and vertical components of the displacement 
vector, 5 denotes a Lagrangian perturbation, and expi(khx — cot) describes the horizontal space 
and time dependence of the perturbations. The horizontal component of the propagation vector, 
kh, is related to the angular degree, £, by 
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where R is the stellar radius. After some manipulation, equations (^) and (||) yield 
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Notice that these two equations are independent of the assumption of adiabaticity. 

The linear, adiabatic wave equation for the fractional, Lagrangian pressure perturbation, 
Sp/p, follows from combining equations (|l|), (||), and (0), 
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The advantages of choosing Sp/p as the dependent variable will become apparent as we proceed. 

The appendix is devoted to deducing the properties of g-modes from the wave equation (||). 
Here we summarize results that are needed in the main body of the paper. A g-mode's cavity 



6 Readers are referred to Unno et al. (1989) for detailed derivations of these equations. 
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coincides with the radial interval within which ui is smaller than both N and Lg. Inside this cavity 
the mode's velocity field is relatively incompressible. The lower boundary of each g-mode's cavity 
is set by the condition to = N. For modes of sufficiently high frequency, the relation u> = is 
satisfied in the radiative interior, and it determines the location of the upper boundary. Otherwise, 
the upper boundary falls in the transition layer between the radiative interior and the convective 
envelope where uj = N . In this paper, we restrict consideration to modes which do not propagate 
immediately below the convection zone since these include the overstable modes. To order of 
magnitude, the condition uj ~ is satisfied at z w ~ to 2 /(gkf t ). If we denote the depth of the first 
radial node of Sp/p as z\, and the depth at the bottom of the surface convection zone as Zb, then 
Z\ ~ z w provided z UJ > z^- 

The fractional Lagrangian pressure perturbation, dp/p, is almost constant in the convection 
zone (z < Zb), and declines smoothly to zero in the upper evanescent region (zb < z < z\). Its 
envelope scales with depth in the cavity so as to maintain a constant vertical energy flux. An 
equal amount of energy is stored between each pair of consecutive radial nodes. For modes having 
z<jj > Zb, the photospheric value of the normalized eigenfunction is given by (5p/p) p \ l ~ l/(nr w L) 1 / 2 . 

3. Perturbations Associated with G-Mode Pulsations 

Here we derive expressions relating various perturbation quantities associated with g-modes.[] 
These are used in §|] to evaluate the driving and damping in various parts of the star. Since 
theoretical details are of limited interest, we summarize selected results at the end of this section. 

We start with the radiative interior and proceed outward through the convection zone to the 
photosphere. A new symbol, A, is introduced to denote variations associated with a g-mode at a 
particular level, such as the photosphere, or within a particular layer, such as the convection zone. 
These variations are not to be confused with either Lagrangian variations which we denote by 8, 
or Eulerian variations. 

This section is replete with thermodynamic derivatives. In most instances we take p and 
T as our independent variables. Unless specified otherwise, it is implicitly assumed that partial 
derivatives with respect to one are taken with the other held constant. We adopt the shorthand 
notation: p p = dlnp/dln p, px = dlnp/dlnT, s p = ds/dlnp, st = ds/dlnT, k p = dlnn/dlnp, 
ht = dlnn/dlnT. However, we use the symbols p s = (din p/ds) p and c p = T(ds/dT) p . These 
thermodynamic quantities evaluated for a model DA white dwarf are displayed in Fig. [2|. 



7 When reading this section, it is important to bear in mind that all perturbation quantities have an implicit time 
and horizontal space dependence of the form expi(khx — ut). 
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Fig. 1. — The squares of the radian Brunt- Vaisala frequency, iV 2 , and Lamb frequency, L 2 , in 
s~ 2 as functions of pressure in dynecm~ 2 for Bradley's white dwarf model with T e g = 12, 000 K. 
iV 2 is negative in the convection zone (logp < 9.3), and drops with increasing depth to near 
zero in the degenerate interior (logp > 23). The little bump at logp ~ 19 is due to the 
sharp compositional transition from hydrogen to helium. The long dashed line is our analytic 
approximation, iV 2 ~ (g/z). Notice that it worsens as degeneracy becomes important. G-modes 
with frequency u propagate in regions where both to < N and u < Lg. These regions are depicted 
by horizontal dashed lines for two I = 1 modes having periods of 134 s and 1200 s respectively. The 
former has its upper cavity terminated where u> = Lg, while the latter propagates until immediately 
below the convection zone where uj = N . 
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Fig. 2. — Values of various dimensionless thermodynamic quantities as a function of pressure in 
dyne cm" 2 for Bradley's white dwarf model with T e g- = 12, 300 K. The convection zone extends 
from the surface down to logp ~ 8.5. Both st and kt change significantly with depth in the 
upper layers where hydrogen is partially ionized. The fully ionized region has mean molecular 
weight fi = 1/2. Other thermodynamic quantities can be obtained from those plotted by using the 
relations: p s = Pt/(pts p - p p s T ), c p = (s T p p - s p p T )/p p , and Ti = c 2 s p/p = (p p s T - p T s p ) / 's T . 
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3.1. Radiative Interior 

The radiative flux in an optically thick region obeys the diffusion equation, 

„ 4a dT A 



3np dz ' 

where a is the Stefan-Boltzmann constant. The Lagrangian perturbation of the flux takes the 
form 

SF 6p .ST d£ z fdlnTy 1 d ( ST\ 



0) 



Next we express 5F/F in terms of 5p/p within the quasiadiabatic approximation. To 
accomplish this, we decompose Sp/p and 5T/T into adiabatic and nonadiabatic components by 
means of the thermodynamic identities 

5p p Sp 3 6p Ss 

— = - 2 + p s 5s^- -, 11) 

p c%p p 5 p 5 

and 

6T _ _ s p P Sp Ss ^2 5p 8s 
1 st c%p p c p op 5 
where the coefficients are evaluated as appropriate for fully ionized hydrogen plasma (cf. Fig. |2|). 



Substitution of the adiabatic parts of these expressions into equation (10) is carried out separately 
for the upper evanescent layer and the g-mode cavity. 



3.1.1. Upper Evanescent Layer 

Within the region z^ < z < z^, 5p/p « —ikh^h, and it varies on spatial scale z u (see the 
appendix). Thus equation (|l]) leads to 

^ 8p_ _8p_ , 13 . 
dz p p 

Moreover, the term 

'dlnTy 1 d /8T\ _ /dlnTy 1 d / s p p 6p\ 
dz J dz \ T J \ dz J dz \st c^p p J ' 

is of order (z/z UJ )(8p/p) and can be discarded]^ Making these approximations, we arrive at 

r p 

where 

A = (3-3«p-2«T) (16) 



8 As shown in Fig. ^| the factors s p ~ —2, st ~ 3, and p/c? s p ~ 3/5 are each nearly constant in the upper part of 
the radiative interior. 
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3.1.2. G-Mode Cavity 

For z > Zu, perturbations vary on the spatial scale of l/k z , where k z « (N/u>)kh > 1/z 
is the vertical component of the wave vector as given by equation ( |A3| ). Here, the dominant 
contributions to 5F/F come from the last two terms in equation (|l0|) . The first transforms to 

— ^ Rd -ik z r ik z Z~, (17) 

dz gp p 

as most easily seen from equation (pi). It then follows that 



6F ik z p dlnp ( dlnT 



dlnT\ 5p ik z N 2 f p \ 2 dlnp 5p 



dlnp / p g \gp J ami p 



F gp dlnT \ dlnp 
where we have set p p ~ 1 and pt ~ 1 as is appropriate for a fully ionized plasma (cf. Fig. |2|). 



3.2. Convective Envelope 

The net absorption of heat by the convective envelope is given by 

AQ = [ dzp—T5s^As b ( dzp—T. (19) 

J cvz J cvz 

Here As;, is the variation of the specific entropy evaluated at the bottom of the convection zone; 
it is neither a Lagrangian nor an Eulerian variation. The chain of argument leading to the second 
expression for AQ requires some discussion. At equilibrium (no pulsation), convection is so 
efficient that ds/dz <^ s/z except in a thin superadiabatic layer just below the photosphere; the 
unperturbed convection zone is nearly isentropic. Moreover, the rapid response of the convective 
motions ensures that dbsjdz <C 5s/ z; the vertical component of the perturbed entropy gradient is 
small during pulsation. Because both ds/dz and d5s/dz w 0, Asb ~ 5s . 

To relate 5s with the flux perturbation, we treat convection by means of a crude mixing-length 
model. At equilibrium, the convective flux is set equal to[] 

& tt pk B T ds 

F ~ aH p v cv — , (20) 

m p dz 

where a is the ratio of mixing length to the pressure scale height. The convective velocity, v cv , 
satisfies^ 

vl ~ -(aH p ) 2 g Ps ^ = - (aH p N) 2 . (21) 

9 Where the convection is efficient, F is comparable to the total flux. 
10 Note that p s =pt/(pts p — p p sr) < 0. 
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Eliminating ds/dz and solving for v cv , we find 

F 

Reversing the procedure and solving for ds/dz yields 
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where s p h is the entropy at the photosphere. The dimensionless factor / ~ a -4 / 3 is of order unity. 
The right hand side of equation (23) is evaluated at the photosphere, since the entropy jump is 
concentrated immediately below it. 

We assume that the mixing length model applies during pulsation, as is consistent with the 
short convective turn-over time. Thus we write the variation of s& — s p h associated with a g-mode 

as 

A(s b - s ph ) _ ^A-Fph 
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(25) 



The dimensionless number C is to be evaluated at the photosphere. In practice, only a few terms 
in the complicated expression for C make significant contributions. In arriving at equations ( |24| ) 
and (^), we implicitly assume that convection carries the entire stellar flux, and that / does not 
vary during the pulsational cycle. 



3.3. Photosphere 

The photospheric temperature and pressure are determined by 



- p\ 1/4 

r«(-) , (26) 



a 



and 



if- 



Provided most of the flux in the photosphere is carried by radiation, the temperature gradient 
is set by the equation of radiative diffusion (eq. f§).[3 Combining equations (||), (p6|), and (| 



n Applymg the diffusion equation at the photosphere is a crude approximation. 
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yields 

On the other hand, the adiabatic gradient, 
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(28) 



(29) 



PTSp ~ PpS T 

using the numbers in Fig. S Comparing the photospheric and adiabatic temperature gradients, 
we find that the convection zone extends up to the photosphere for DA white dwarfs which lie 
inside the instability strip. This is confirmed by comparison with numerical models (e.g., Bradley 
1996j Wu |i~997D . However, b ecause the photospheric density is low, convection is inefficient in this 



region, and radiation carries the bulk of the flux. 

Next we relate changes in the thermodynamic variables at the photosphere to changes in the 
emergent flux. Thus 

AT _ 1 AF ph 
T ~ 4 F 

is an immediate consequence of equation (f26|). Expressions for Ap/ p and Ap/p follow in a 
straightforward manner from the thermodynamic identity, 
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and from the perturbed form of equation (^7|), 
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It is then a simple step to show that 
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3.4. Putting It All Together 

We begin by collecting a few key equations obtained in previous sections. These are then 
combined to determine how the entropy and flux variations in the convection zone and photosphere 
depend upon the pressure perturbation associated with the g-mode. 



3.4-1- Key Equations 

The flux perturbation entering the bottom of the convective envelope from the radiative 
interior satisfies (eq. |15tl ), 

(37, 

The photospheric entropy varies with the flux perturbation emerging from the convection zone as 
(eq. @]) 

As ph ^B^. (38) 

The variation in the superadiabatic entropy jump is related to this photospheric flux perturbation 
by (eq. ||) 

A(s b -s ph )^C^. (39) 

The dimensionless numbers A, B and C are all positive; A is evaluated at the boundary between 
the radiative interior and the convective envelope, whereas B and C are computed at the 
photosphere. 

The net heat variation, AQ, and specific entropy variation, As;,, of the convective envelope 
are connected by equation (19), 

AQ « F U As b , (40) 

where the thermal time constant 

n^-! dz^* p -^. (41) 

F J cvz m p 7F 
The set of key equation is completed by the relation 

^ = AF b - AF ph . (42) 
at 

We ignore horizontal heat transport for good reasons: transport by radiation is completely 
negligible because k^z^ <g 1; turbulent diffusion acts to diminish Asj, but only at the tiny rate 
k\zv cw <C uj. 

For compactness of notation, we define a new thermal time constant, r c , by 

t c = (B + C)n- (43) 

The physical relation between r c and r& is discussed in §|B]. Here we merely note that r c is generally 
an order of magnitude or more larger than Tb- 
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3.4-2. Implications Of Key Equations 



Taken together, the five homogeneous equations (p7\), (pq), (p9[), (|40|), and ([42]) enable us 
to express the five quantities As&, As p h, AFb/F, AF p h/F, and AQ in terms of (5p/p)&.P| The 
principal results from this section read: 

^„4Z+£i(*), (44, 



1 - iujT c \p J b 
A, ph M (*) , (45) 

1 - IUJT C \ p J b 

AF >-A( S A, (46) 



F ' VP 



6 



^ « _ ^— (47) 

Ag«^(*). (48) 

1 - ZWTc V P J b 



4. Driving and Damping 

The produce of this section is the proof that a g-mode is linearly overstable in a ZZ Ceti star 
provided the convective envelope is thick enough, or more precisely, provided ujt c ^ l.fj 

The time-averaged rate of change in a mode's energy, 

uj f 7 1 dE 

is obtained from the work integral. Useful forms for 7 areP| 

7,2 / , f R , k B rrT ,d6s co T f , f R , 5T d fSF\ , , 

Regions of driving and damping are associated with positive and negative values of the above 
integrand. The evaluation of 7 for modes with z^ > z b is simplified by the near constancy of 5p/p 
for z <C z u . 



12 We include AF^/F in this list for completeness although it is already expressed in terms of (5p/p)b by equation 



1 The qualification that > zj, is necessary here. 

14 These assume that the eigenfunction is normalized such that E = 1. 
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4.1. Radiative Damping 

In the upper evanescent region of the radiative interior, the perturbed flux varies in phase 
with the pressure perturbation because adiabatic compression causes the opacity to decrease.^]. 
Since 5p/p declines from close to its surface value at Zj, to zero near , and then oscillates 
with rapidly declining amplitude at greater depth, the upper evanescent layer loses heat during 
compression and thus contributes most to mode damping. 

To evaluate the radiative damping rate, it proves convenient to use the second form for 7 
given by equation (|50|). We consider the contributions to 7 ra a from the upper evanescent layer, 
7 ru , and from the propagating cavity, 7 n , separately. 

The contribution from Zb < z < z u , obtained with the aid of equation (|l~^) for 5T/T and 
equation fl46[ ) for 5F/F, reads 

fAL\ /5p\ 2 L f z " , /5p\ 2 dA . . 

-T7T - / ^ - -H- (51) 



10 y V v J b 10 7 Zb vp/ dz 

The first term on the right hand side generally dominates over the second one as dp/p declines 
with depth and A does not vary significantly in this region. 



For z > z^, we again use equation ( |12| ) for 5T/T, but now substitute equation (|lg|) for AF/F 
to arrive at 

125ff 3 a; 2 y Za; dlnTVp 
Appeal to equation (|A8|) giving the WKB envelope relation for <5p/p establishes that the integrand 
peaks close to Thus, 



7 „ / dziV 4 c 4 ^ ^ . (52) 



— (f):(?):~-K?):- 

where in this context (Sp/p)^ represents the magnitude of the WKB envelope evaluated at z u . 
Since this magnitude is significantly smaller than (Sp/p)b, the contribution to 7 from the g-mode 
cavity is negligible provided > Zb- 

To a fair approximation, 7 ra 4 may be set equal to the first term on the right-hand side of 
equation (^lj). Thus 



4.2. Convective Driving 

The perturbed flux which exits the radiative interior enters the bottom of the convection zone. 
There it is almost instantaneously distributed so as to maintain the vertical entropy gradient near 



15 A > in equation @ 
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zero. Because the convection zone gains heat during compression, it is the seat of mode driving. 

To evaluate the rate of convective driving we substitute equations (12) and (f44|) into the first 
form for 7 given by equation (50). It is apparent that the net contribution comes entirely from the 
adiabatic part of 5T/T (see eq. [12|). Since the integrand is strongly weighted toward the bottom 
of the convection zone, we evaluate all quantities there and arrive at 



(J 7V, 



1 + K) 2 



AL 



(55) 



Since A is positive, so is 7 C 



4.3. Turbulent Damping 

Damping due to turbulent viscosity acting on the velocity shear in the convection zone as well 
as the convective overshoot region is estimated as 

7visc « -4?ri2V £ dzpu . (56) 

Here, v ~ v cv H p is the turbulent viscosity, and —iuod^/dz is the dominant component of velocity 
gradient. 

Adiabatic perturbations in an isentropic region are irrotational. Consequently, the velocity 
shear in an isentropic convection zone would be very small and lead to negligible turbulent 
damping. However, the mean entropy gradient and the horizontal gradient of the entropy 
perturbation represent significant departures from isentropy. In a future paper we demonstrate 
that turbulent viscosity suppresses the production of velocity gradients in the convection zone to 
the extent that its effect on damping is negligible in comparison to radiative damping. However, 
turbulent damping in the region of convective overshoot at the top of the radiative interior 
stabilizes long period modes near the red edge of the ZZ Ceti instability strip. 



4.4. Net Driving 



The net driving rate follows from combining equations (54) and (| 

(wr c ) 2 - 1 (AL\ f5p^ 2 



7net 



(57) 



K) 2 + 1 V 10 j \ P/b 

Driving exceeds damping by a factor of two in the limit ujt c S> 1- Substitution of the normalization 



relation given by equation (A28) for 5p/p yields the more revealing form 



1 



7net 



K) 2 + 1 



1 



(58) 
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Overstability occurs if lot c > Following Brickhill, we refer to this excitation mechanism as 
convective driving. 



5. Discussion 
5.1. Time-Scales 

Three time-scales are relevant for convective driving in DAVs. The first is the period of an 
overstable g-mode, P = 2it/uj, which is typically of order a few hundred seconds. The second is the 
dynamical time constant, t cv ~ H p /v cv , on which convective motions respond to perturbations; 
t C v < 1 s throughout the convection zones of even the coolest ZZ Cetis. This is why the convective 
motions adjust to the instantaneous pulsational state. The third is the thermal time constant, r c , 
during which the convection zone can bottle up flux perturbations that enter it from below. 

Given the central role of r c , we elaborate on its relation both to t cv and to the more 
conventional definition of thermal time constant, r t h, at depth z. The latter is the heat capacity 
of the material above that depth divided by the luminosity. In a plane parallel, fully ionized 
atmosphere this is equivalent to 



_ 1 f z , pk B _ 5pz 

7th = -=; / UZ C D — PS — — . (59) 

F Jo P m n T 7F 



Appeal to equation (^) establishes that inside the convection zone 



— «1. (60) 

Tth \C S J 

Now r c = (B + C)^, where is defined by equation (|4l|). To the extent that c p ~ 5 is constant in 
the convection zone, ph r t h/5, where the latter is evaluated at z\,. 

Next we address the relation between r c and r&. Here we are concerned with the relatively 
large value of B + C, typically about 20 for DAVs-P 7 ! Recall from equations ( f44|) and ( ^7|) that 



B + C 



(61) 



So the photosphere and superadiabatic layer add an insulating blanket on top of the convection 
zone. The large value of B follows because the photospheres of DAVs are composed of lightly 
ionized hydrogen. In this state, the values of kt and st are both large and positive; typical values 
in the middle of the instability strip are ki>«6 and st ~ 24. The large and positive kt arises 
because the population of hydrogen atoms in excited states which the ambient radiation field 



16 In the quasiadiabatic limit for modes with > z\,. 
17 In our models, B and C have comparable value. 
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can photoionize increases exponentially with increasing T. The large and positive st occurs as 
the ionization fraction increases exponentially with increasing T, and the much larger entropy 
contributed by a free as compared to a bound electron. The large and positive value of C reflects 
the increase in entropy gradient that accompanies an increase in convective flux. It is obtained 
from mixing length theory with an unperturbed mixing length. 



5.2. Validity of the Quasiadiabatic Approximation 

The validity of the quasiadiabatic approximation requires that the nonadiabatic parts of the 
expressions for 5p/p and ST/T, as given by equations ( |Tl| ) and (12), be small in comparison to the 
adiabatic parts. Thus the ratio 

n - s Wp (62) 

is a quantitative measure of nonadiabaticity. We estimate TZ na , for the radiative interior and the 
convection zone. 

We calculate 5s in the radiative interior from 

iF m p d {5F\ 

6s ^-^fT z {-¥)- (63) 

In the upper evanescent layer, z\, < z < z u , this leads to 



l^nal ~ — , (64) 



whereas in the propagating cavity, z > ZoJ , we find 



\n na \ ~ — ( — ) . (65) 



Nonadiabatic effects in the radiative interior are maximal at z = z^, where 



l^nal ~ — , (66) 



since t^/z increases with depth. 



The measure of nonadiabaticity in the convection zone is given by equation fl44|). Since 
uJTf, > 1 is required for the validity of the quasiadiabatic approximation in the radiative zone, we 
restrict consideration to the limiting case lot c 3> 1. In this limit, equation yields 



|ftna| ~ — , (67) 
UJT b 



which is identical to the value arrived at for the radiative zone in equation (I 
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The requirement ujTb ;> 1 for the validity of the quasiadiabatic approximation severely limits 
the applicability of the current investigation. The perturbed flux at the photosphere is related to 
that at the bottom of the convection zone by 



AF ph 1 AF b 



1 — iujr c F 



(68) 



Since r c is at least an order of magnitude larger than r^, modes with ujt^ ;> 1 are likely to exhibit 
small photometric variations. However, this may not render them undetectable because their 
horizontal velocity perturbations pass undiminished through the convection zone. 



5.3. Brickhill's Papers 



Our investigation is closely related to studies of ZZ Cetis by Brickhill ( [L983|Jl990| , 1991a ) 



Brickhill recognized that the convective flux must respond to the instantaneous pulsational state. 
To determine the manner in which the convection zone changes during a pulsational cycle, he 
compared equilibrium stellar models covering a narrow range of effective temperature. Brickhill 
provided a physical description of convective driving and obtained an overstability criterion 
equivalent to ours.[^] Moreover, he recognized that the convection zone reduces the perturbed flux 
and delays its phase. 

Our excuses for revisiting this topic are that Brickhill's papers are not widely appreciated, 
that our approach is different from his, and that our paper provides the foundation for future 
papers which will examine issues beyond those he treated. 



A. Scaling Relations for Gravity-Modes 

The appendix contains derivations of a number of scaling relations appropriate to g-modes 
in ZZ Cetis. These relations are applied in the main text and in our subsequent papers on white 
dwarf pulsations. We obtain approximate expressions for dispersion relations, eigenfunctions, and 
normalization constants. 



A.l. Properties of Gravity Waves 

Substituting the WKB ansatz 



— = *4,exp ( / I (It /,•- j , ( Al i 



3 Our time constant r c is equivalent to the quantity D which Brickhill defined in equation (9) of his 1983 paper. 
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into the wave equation (||), and retaining terms of leading order in k z H p S> 1, we obtain the 
dispersion relation 

(N 2 - uj 2 ) (k 2 c 2 - uj 2 ) 



C 2 UJ 2 



(A2) 



where khC s is the plane parallel equivalent of the Lamb frequency, Lg. 



Gravity waves propagate in regions where uj is smaller than both N and Lg. Far from turning 
points, their WKB dispersion relation simplifies to 



k z ~ i kh- 
uj 

The vertical component of their group velocity is given by 



duj 
dk z 



UJ _|_ 



(A3) 



(A4) 



Terms of next highest order in k z H p yield the WKB amplitude relation for gravity waves; 



A 2 cx 



P oc P 



k z p 2 Np 2 



(A5) 



The amplitude relation expresses the conservation of the vertical flux of wave action. This may 
be confirmed directly by means of an alternate derivation. A straightforward manipulation of the 
linear perturbation equations yields the quadratic conservation law 









at 


5 


at 



+ 



(p 



t\2 



2pc 



2 ~ 9 Iy 



X7 I > d * 



(A6) 



where p' = bp — gp^ z is the Eulerian pressure perturbation. We identify F = p'(d£dt) as the 
quadratic energy flux. The time-averaged magnitude of the vertical component of F, computed to 
lowest order in {k z z)~ x with the aid of equations ([?]) and ([A3]), is found to be 



^ uj 2 Np 2 (Spy 



c/» gp \ p / 



(A7) 



which is in accord with the WKB amplitude relation. Yet another way to look at the energy 
flux is as the energy density transported by the group velocity. The energy density may be 
approximated as oj 2 p&, since kinetic energy contributes half the time-averaged total energy 



density, and |^| S> Application of equations (ph, (|A3j), and (A4) yields 



T uj 2 p£ h v gz 



uj 2 Np 2 /5p^ 2 

gkh gp \p 



(A8) 



which reproduces the result for T given by equation (A7). 

In a spherical star, the conserved quantity is the WKB luminosity, C = Airr 2 ^. 
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A.2. Properties of G-Modes 

A. 2.1. Inside the Propagating Cavity 

The Lamb frequency, Li, decreases with increasing r, whereas the Brunt- Vaisala frequency, 
N, increases until it drops abruptly just below the convection zone and becomes imaginary 
within it (cf. Fig. 1). Compositional transitions may be exceptions to these trends. The WKB 
approximation is violated if these transitions are sharp on the scale of the radial wavelength. When 
these are located in regions of propagation, they are best viewed as separating linked cavities. 

We divide a star's g-modes into two classes, high frequency modes and low frequency 
modes. The former have cavities bounded from above at z w > Zb where to = Li ~ khC s . The 
latter propagate until just below the convection zone at z ~ z^ where u = N. The floor of the 
propagating cavity is set by oj = N for all modes. In this paper we are exclusively concerned with 
high frequency modes. Taking c s ~ (gz) 1 / 2 , we find 



j 1 



Moreover, N ~ (g/z) 1 / 2 ~ c s /z implies k z z^ ~ 1 at z = z u . Thus z^ is to be identified with zi, 
the uppermost node of 5p/p. More generally, 

K ~ — L^, (A10) 

for z > Zuj. 

Expressions for ^ and £ z are given by equations @ and (^). To leading order in k z H p , these 
reduce to 

k z P ( $P\ 



kh gp V p J 



and 

p ( 5p 



gp \p 

Where N/uj ^> 1, we have k z /kh ~S> 1, which implies \£,h/£,z\ 3> 1- In these regions g-modes 
are characterized by nearly horizontal motions. Equations ( |A11| ) and ( A12| ) imply that 



(A12) 



5p/p = —ikhCh — d^z/dz vanishes to leading order in k z H p ^> 1. Thus g-modes are relatively 
incompressible within their cavities. A more precise estimate, obtained from equation ( All| ), is 
Sp/p ~ ikhihiz^/z) 1 / 2 . 



5 10 15 20 

log p 



Fig. 3. — Structure of a mode with n = 9, £ = 1, and a period of 502 s, as a function of pressure 
in dynecm" 2 . The top panel, which is similar to Fig. [j], illustrates how the mode cavity, depicted 
by the dashed line, is formed. For this mode z u > z^. The WKB luminosity measured in erg/ s is 
plotted in the second panel. It is constant inside the cavity and decays outside it. The compositional 
transition between the hydrogen and helium layers has minimal effect on this mode. The lower 
two panels display the depth dependences for the dimensionless 5p/p and ^ measured in cm. The 
numerical values in this figure come from setting ^ = R at the photosphere. 
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A.2.2. Upper Evanescent Layer 



For the purpose of this discussion, we pretend that both p — > and p — * as z — ► 0. Then 
z = is a singular point of the wave equation (H). Since the physical solution is regular at z = 0, 



In the convection zone, where N 



d_ 

dz 



d 

dz 



In 



<5p 



k\z. 



In 



Vcv_ 



and at the top of the radiative interior, where iV 2 ~ (c s /z) 2 , 



d_ 

dz 



In 



<5p 



(A13) 



(A14) 



(A15) 



We see that Sp/p is nearly constant for z C ^ in the upper evanescent layer .P| 

Taking into account the near constancy of Sp/p, equations (||) and (^) imply that to leading 
order in z/z^, 

(A16) 
(A17) 



and 



u; 2 / 5p 
5^1 V V 

Thus both components of the displacement vector are nearly constant for z <C The behavior 
of is more subtle; in principle, it could vary on scale z should p/{c 2 s p) do so. However, in 
practice, p/(c 2 p) exhibits only mild depth variations. Thus, equation (|A16f) yields 



Sp 
P 



-ikh£,h- 



(A18) 



Equation (A18) shows that the relative incompressibility that characterizes propagating g-modes 
does not extend to their evanescent tails. 



A. 3. Global Dispersion Relation 

The global dispersion relation is obtained from 



i t 77 ■■ I dz k z , 



(A19) 



9 That is why we chose to use it as the dependent variable in the wave equation. 
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where z\ is the lower boundary of the mode cavity. Carrying out the integration using N 2 ~ g/z, 
Z\ ~ z u , and equation (|A3|), we obtain u> 2 ~ gkh/n. This relation is a good approximation to the 
g-mode dispersion relation for a polytropic atmosphere. The proportionality 

lo 2 oc ^ (A20) 
n 

provides a satisfactory fit to the high frequency modes in DAV white dwarfs. However, low 
frequency modes penetrate deeply into the interior where the approximation N ~ {g/z) 1 / 2 fails 
because of electron degeneracy (cf. Fig. [I]). As a result of the steep drop in N, Z\ is nearly 
independent of uj for uj < 1CP 2 s _1 . This steepens the dependence of uj on n such that 

woe—. (A21) 
n 

Numerical results for the dispersion relations are shown in Fig. ||. 



A. 4. Normalization of Eigenfunctions 

We conform to standard practice and set 

^j*drr 2 P (e h + ii)=l. (A22) 

The conservation of the WKB energy flux implies that the region between each pair of neighboring 
radial nodes contributes an equal amount to the above integral. 

To achieve a simple analytic result, we take advantage of the following: |^| 3> |£ z | and 
ih ~ Chlph f° r z <^ z w , and the envelope of pv gz £, 2 ~ constant for z ^> z u (cf. eq. (|A8|| ). This 
enables us to write 

Y^U^L/^-l, (A23) 

z J u gz 

where the subscript uj stands for quantities evaluated at z w . Using 

R dz Tin 

— « — , (A24) 

Vgz U 



we arrive at 



where F is the stellar flux, and 



is the thermal time scale at z u 



pz\ 
F , 



(A26) 
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Fig. 4. — Radian frequencies of gravity-modes in s _1 as functions of radial order, n, and spherical 
degree, I. These eigenvalues are computed using Bradley's white dwarf model with T e g = 12, 000 K. 
The upper panel is for n = 1 modes with various angular degrees, while the lower panel is for I = 1 
modes of different radial orders. The global dispersion relations given by equations (|A20|) and 



(A21) fit well for modes of low and high order respectively. We adopt these empirical laws in our 
analytical studies. 
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Putting these relations together, we obtain 



7,2 c 2 1 



1 



where L = AirR F is the stellar luminosity. Equation (A16) then implies 



5_p\ 2 
P / P h 



1 



This is an appropriate normalization formula for modes having z w > zj,. 



(A27) 



(A28) 
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Fig. 5. — Normalized surface amplitudes as well as driving/damping rates for g-modes in white 
dwarfs. We compare numerical values obtained using Bradley's model with T e g- = 12, 000 K 
{jb = 300 s) with analytic estimates. Both properties depend much more strongly on n than 
on l\ hence the choice of n as the abscissa. The upper panel plots as triangles numerical values of 
(8p/p)ph in gm -1 cm~ 2 s 2 for modes with £ = 1 to 10. The eigenfunctions are normalized according 
to equation ( A22j ). The solid line represents the analytical estimate from equation ( A28| ) for modes 
with z u > Zb, where is chosen to be 2>ituj 2 / (T.gk'fy. The lower panel displays numerical radiative 
damping rates as triangles, together with the analytic estimate from equation (|54|) as a solid line. 



